Simulation of transient energy distributions in sub-ns streamer formation

Breakdown and streamer formation is simulated in atmospheric pressure nitrogen for a 2D planar electrode system. A PIC code with multigrid potential solver is used to simulate the evolution of the non-equilibrium ionization front on sub-nanosecond timescales. The ion and electron energy distributions are computed, accounting for the inclusion of inelastic scattering of electrons, and collisionally excited metastable production and ionization. Of particular interest is the increased production of metastable and low-energy ions and electrons when the applied field is reversed during the progress of the ionization front, giving insight into the improved species yields in nanosecond pulsed systems

Numerical modelling of low temperature plasma

The intention of this thesis is to gain a better understanding of basic physical processes occurring in low temperature plasmas. This is achieved by applying both analytic and numerical models. Low temperature plasmas are found in both technological and astrophysical contexts. Three different situations are investigated: an instability in electronegative plasmas; electron avalanches during plasma initiation; and a phenomenon called the Critical Ionisation Velocity interaction. Industrial plasma discharges with electronegative gases are found to be unstable in certain conditions. Fluctuations in light emission, particle number densities and potential are observed. The instability has been reproduced in a variety of experiments. Reports from the experiments are discussed to characterise the key features of the instability. An, as yet un-considered, physical process that could explain the instability is introduced. The instability relies on the plasma's transparency to the electric field. This mechanism is investigated using simple zero-dimensional numerical and analytic models. The results from the models are compared to experimental results. The calculated frequencies are in good agreement with the experimental measurements. This shows that the instability mechanism described here is relevant. For the remaining two problems a three-dimensional particle model is constructed. This model calculates the trajectories of each individual particle. The potential field is solved self-consistently on a computational mesh. Poisson's equation is solved using a Multigrid technique. This iterative solution method uses many grids, of different resolutions, to smooth the error on all spatial scales. The mathematical foundation and details of the components of the Multigrid method are presented. Several test cases where analytic solutions of Poisson's equation exist are used to determine the accuracy of the solver. The implemented solver is found to be both efficient and accurate. Collisions are vitally important to the evolution of plasmas. The chemistry resulting from collisions is the reason why plasmas are so useful in technological applications. Electron collisions are included in the particle model using a Monte-Carlo technique. A basic method is given and several improvements are described. The most efficient combination of improvements is determined through a series of test cases. The error resulting from the collision selection process is characterised. Technological plasmas are formed from the electrical breakdown of a neutral gas. At atmospheric pressure the breakdown occurs as an electron avalanche. The particle model is used to simulate the nanosecond evolution of the avalanche from a single electron-ion pair. Special attention is paid to the inelastic collisions and the creation of metastables. The inelastic losses are used to estimate the photon emission from the electron avalanche. The Critical Ionisation Velocity phenomena is investigated using the particle model. When a neutral gas streams across a magnetised plasma the ionisation rate increases rapidly if the speed of the neutrals exceeds a critical value. Collisions between neutrals and positive ions create pockets of unbalanced negative charge. Electrons in these pockets are accelerated by their potential field and can reach energies capable of ionisation. The evolution of such an electron overdensity is simulated and their energy gain under different density and magnetic field conditions is calculated. The results from the simulation may explain the discrepancy between laboratory and space experiments

Strangeness contribution to the vector and axial form factors of the nucleon

The strangeness contribution to the vector and axial form factors of the nucleon is presented for momentum transfers in the range $0.45<Q^2<1.0$ GeV$^2$. The results are obtained via a combined analysis of forward-scattering parity-violating elastic $\vec{e}p$ asymmetry data from the $G^0$ and HAPPEx experiments at Jefferson Lab, and elastic $\nu p$ and $\bar{\nu} p$ scattering data from Experiment 734 at Brookhaven National Laboratory. The parity-violating asymmetries measured in elastic $\vec{e}p$ scattering at forward angles establish a relationship between the strange vector form factors $G_E^s$ and $G_M^s$, with little sensitivity to the strange axial form factor $G_A^s$. On the other hand, elastic neutrino scattering at low $Q^2$ is dominated by the axial form factor, with still some significant sensitivity to the vector form factors as well. The combination of the two data sets allows the simultaneous extraction of $G_E^s$, $G_M^s$, and $G_A^s$ over a significant range of $Q^2$ for the very first time.Comment: 3 pages, 1 figure, will appear in AIP Conference Proceedings for PANIC 200

An Energy-Minimization Finite-Element Approach for the Frank-Oseen Model of Nematic Liquid Crystals: Continuum and Discrete Analysis

This paper outlines an energy-minimization finite-element approach to the computational modeling of equilibrium configurations for nematic liquid crystals under free elastic effects. The method targets minimization of the system free energy based on the Frank-Oseen free-energy model. Solutions to the intermediate discretized free elastic linearizations are shown to exist generally and are unique under certain assumptions. This requires proving continuity, coercivity, and weak coercivity for the accompanying appropriate bilinear forms within a mixed finite-element framework. Error analysis demonstrates that the method constitutes a convergent scheme. Numerical experiments are performed for problems with a range of physical parameters as well as simple and patterned boundary conditions. The resulting algorithm accurately handles heterogeneous constant coefficients and effectively resolves configurations resulting from complicated boundary conditions relevant in ongoing research.Comment: 31 pages, 3 figures, 3 table

X - Ray Flares and Their Connection With Prompt Emission in GRBs

We use a wavelet technique to investigate the time variations in the light curves from a sample of GRBs detected by Fermi and Swift. We focus primarily on the behavior of the flaring region of Swift-XRT light curves in order to explore connections between variability time scales and pulse parameters (such as rise and decay times, widths, strengths, and separation distributions) and spectral lags. Tight correlations between some of these temporal features suggest a common origin for the production of X-ray flares and the prompt emission.Comment: 7th Huntsville Gamma-Ray Burst Symposium, GRB 2013: paper 15 in eConf Proceedings C130414